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ABSTRACT 


This  paper  considers  the  problem  of  maximizing  the  energy  or  average 

power  transfer  from  a  nonlinear  dynamic  source.  The  main 

theorem  includes  as  special  cases  the  standard  linear  result 

=  Y*  ^  and  a  recent  finding  for  nonlinear  resistive  net- 
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works.  An  operator  equation  for  the  optimal  output  voltage  v(*)  is 
derived,  and  a  numerical  method  for  solving  it  is  given. 
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This  paper  addresses  the  problem  of  extracting  the  maximum  energy  or 
average  power  from  a  nonlinear  dynamic  source  with  the  topology  shown  in 
Fig.  1.  The  main  results  are  i)  a  simple  iterative  scheme  for  finding  the 
optimal  output  voltage  waveform  v(*)  for  any  given  current  source  waveform 
i^(*);  and  ii)  an  expression  for  the  (noncausal )  optimal  load  admittance 
operator  in  terms  of  the  source  admittance.  The  first  result  can  be  useful 
in  engineering  practice  because  it  specifies  the  optimal  performance  that 
is  possible  in  principle  and  because  v(*)  itself  is  a  concrete  design  goal: 
any  load  for  which  the  output  voltage  closely  approximates  v(*)  will  absorb 
nearly  the  maximum  possible  energy  or  average  power.  The  second  result  has 
no  immediate  impact  on  applications  because  the  optimal  load  is  noncausal.  But 
it  has  some  theoretical  significance  because  it  general  izes  and  unifies  previous  work: 
the  1-port  versions  of  thestandard  result  ^source^®**  ‘^y^^mic  linear^ 

time-invariant  (LTI)  systems  [3]  and  a  recent  theorem  for  resistive  nonlinear 
systems  [4]  fall  out  as  special  cases. 

The  body  of  this  paper  addresses  the  topic  rigorously  by  giving  sufficient 
conditions  for  existence,  uniqueness  and  global  optimality  of  the  network 
solution,  along  with  a  convergence  proof  for  the  iterative  algorithm.  But 
the  remainder  of  this  introduction  is  utterly  nonrigorous  and  enables  the 
reader  to  sample  the  results  before  (or  instead  of)  delving  into  the  mathe¬ 
matics  used  to  establish  them. 

1.1)  Informal  Description  of  Results 

Restricting  attention  to  sources  with  the  special  "Norton  form"  topology 


1.  Reference  [3]  actually  deals  with  the  dual  network,  where  the  source  appears 
in  Thevenin  form. 


-2- 


shown  in  Fig.  1  makes  it  possible  to  derive  a  load  admittance  operator 

that  is  optimal  for  any  current  source  waveform  i5(’)-  The  description  of 

Gopt  involves  the  linearized  behavior  of  the  source  admittance  operator  F 

about  any  nominal  voltage  waveform  v„^„(*).  i.e., 

nom  ’ 


"(''.on,'-*) 


(t,T) 6v  (t)  dr 


(1) 


nom 


to  first  order  in  5v(*).  where  the  second  term  on  the  right  is  linear  in  5v(*) 
and,  in  general,  time-varying.  Under  certain  assumptions  on  F  described 
later,  the  optimal  load  turns  out  to  be 


00 

Gopt(v('))  =  I  h^(T,t) V  (x) dx 


(2) 


Note  that: 

i)  in  the  causal  LTI  case  where  F  is  characterized  by  an  impulse  response 
hy(t,T)  *  h(t-T)  with  Fourier  transform  ,  Eq.  (2)  describes 

an  anti-causal  LTI  load  with  impulse  response  h(x-t)  and  thus  an  admit- 
tance  ,  in  agreement  with  the  classical  result  [3]; 

ii)  if  the  source  is  a  nonlinear  resistor,  i  =  f(v),  then  the  optimal  load 
is  a  nonlinear  resistor  i  =  v*f'(v),  in  agreement  with  [4]; 

iii)  represents  the  behavior  of  F  linearized  about  the  independent 

variable  v(*)  as  a  nominal  input,  so  G^p^  is  nonlinear  in  general; 

iv)  for  causal  F,  G^p^  is  unrealizable  in  all  but  the  purely  resistive  case 
since  the  roles  of  t  and  x  in  Eq.  (2)  are  reversed  from  Eq.  (1); 

v)  Eq.  (2)  is  exact  for  large-signal  behavior,  despite  the  fact  that  Eq. 

(1)  is  only  a  valid  approximation  for  small  6v(*)  ; 

vi)  first  appearances  notwithstanding,  G^^p^  is  time-invariant  if  F  is;  and 
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vii)  the  optimal  load  admittance  operator  is  a  linear  function  of  the  source 
admittance.  Thus  the  optimal  load  for  a  parallel  connection  of  source 
admittances  is  the  parallel  connection  of  the  optimal  loads  for  each 
source  separately. 


Given  The  optimal  voltage  v(*)  is  the  solution  to  the  network 

in  Fig.  1 ,  i . e • , 

i5(0  =  F(W-))  +  ^ 

where  H  is  the  combined  admittance  loading  the  current  source.  And  Eq. 

(3)  can  be  solved  in  practice  by  the  iterative  procedure, 

Vj^.i(*)  =  Vj.{-)  +  r(i^(-)  -  H(Vj.(-))  ,  (4) 


which  is  guaranteed  to  converge  given  any  initial  guess  Vq(*)  if  r  >  0 
is  sufficiently  small  and  certain  technical  conditions  are  satisfied. 
The  solution  v(*)  gives  both  the  circuit  behavior  when  optimally  loaded 
and  the  performance  bound  for  average  power  transfer,  which  can  be  cal¬ 
culated  as 


max 


I 

T 


v(t)[G^p^(v(-))](t)dt 


(5) 


for  a  drive  with  period  T. 

In  general,  no  causal  load  will  produce  the  optimal  output  voltage 
for  every  i5(').  but  in  practice  one  frequently  encounters  the  more 
restricted  problem  of  maximizing  power  transfer  for  some  single 
In  this  case  there  may  be  a  variety  of  causal  loads  that  will  do  the  job, 
i.e.,  load  admittances  such  that  Eq.  (3)  and 


i^C-)  =  F(v(.)) 

have  the  same  solution  v(*)  for  the  particular  ijl*)  of  interest.  In 


(6) 


designing  a  circuit  to  realize  »  ''(•)  can  serve  as  a  design  goal 

and  as  a  performance  standard. 

fnaX 

1.2)  Example 

Suppose  the  source  takes  the  specific  form  shown  in  Fig.  2,  with 

k-1 

the  resistor  curves  given  by  i  =  “  vjvl  ,  k  =  1,2,3,  as  shown  in 

Fig.  3.  Then: 

00 

[Fk(v(-)):(t,  •  g|^{v(t))  +  I  »(T)dT  ,  (7) 

•  00 

where  l(x)  is  the  unit  step  function  vanishing  for  x  <  0.  The  optimal 
load  admittance  is 


[GQp^_(v(-))](t)  =  v(t)-g'(v(t))  +  [  l(T-t)e’^^-^^(T)dT 


(8) 


One  can  check  that  the  resistor  curves  are,  in  fact,  continuously  dif¬ 
ferentiable,  with  derivatives  given  by 

gk(v)  =  klvl"’^  ,  k  =  1,2,3  .  (9) 

The  optimal  output  voltage  v(*)  was  determined  by  numerically  carrying 
out  the  iterative  procedure  (Eq.  4),  which  in  this  case  takes  the 
form 

J  *  CO 


6  sin(t)  -  (k+l)Vj(t) |Vj(t) 


+  Vj(t),  k  =  1,2,3 


,k-l 


e~'^"^' v.(T)dT 

J 


(10) 


Since  represents  a  linear  resistor,  it  follows  from  the  tradi¬ 
tional  linear  theorem  that  v{t)  =  2 sin  (t)  for  k=l,  in  agreement  with 
the  numerical  solution.  Note  that  the  instantaneous  current  drawn  by 
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the  nonlinear  source  resistor  increases  in  magnitude  with  k  for  |v|  >  1 
but  decreases  for  !v|  <  1.  Thus  it  is  intuitively  reasonable  that  with 
increasing  k  the  optimal  output  spends  a  progressively  greater  percen¬ 
tage  of  time  in  the  region  |vl  <  1,  as  seen  in  Fig.  4. 


1.3)  Generality  and  Limitations 


The  results  in  this  paper  extend  in  a  straightforward  way  to  the  dual 
case,  i.e.,  "Th^venin  form"  topology,  consisting  of  an  independent 
voltage  source  in  series  with  a  circuit  element  characterized  by  an 
impedance  operator.  The  extension  to  multiport  sources  is  also  straight 
forward. 

For  source  networks  with  the  topology  shown,  the  key  restriction 
is  that  H  be  monotone  increasing,  i.e.,  incrementally  passive.  But  the 
assumed  topology  is  perhaps  a  greater  restriction,  since  the  results 
have  been  shown  to  depend  on  this  topology  in  a  fundamental  way  [5] 
and  since  there  is  no  general  nonlinear  analog  of  the  transformation 
used  to  put  any  linear  circuit  in  Thevenin  or  Norton  form. 

The  difficult  nonlinear  version  of  the  linear  "broadband  match¬ 
ing  problem"  [6-9],  in  which  the  goal  is  to  choose  a  causal  load  ad¬ 
mittance  that  optimizes  power  transfer  for  a  set  of  inputs,  is  not 
addressed  in  this  paper  and  remains  entirely  open,  to  the  best  of  the 
author's  knowledge. 


2.1)  Notation  and  Definitions 

Let  L  be  any  real  inner  product  space  and  L  any  linear  subspace  of  L. 

An  operator  F:  L  -*■  L  is  said  to  be 

a)  strictly  increasing  if 

<'F(y)  -  F(x),  y-x)  >  0,  Vx  ^  y  e  L  ,  (11) 

b)  uniformly  increasing  if  for  some  6  >  0  , 

{F(y)  -  F(x),  y-x^  >  6  ||y-x|!  Vx.y  e  C  ,  (12) 

c)  Lipschitz  continuous  if  for  some  K  1  0  , 

!|F(y)  -  F(x)(|  1  K||y-xl|,  Vx.ye  L  .  (13) 

Let  L,  L'  be  any  real  inner  product  spaces,  and  L(L,L')  denote  the  space 

of  continuous  linear  maps  from  L  to  L',  with  the  operator  norm  [10,  p.316]. 

For  AeL(L,L'),  let  A^^'^  denote  the  adjoint  of  A. 

2 

The  Hilbert  space  L  is  the  set  of  all  measurable,  square-integrable 
functions  x:  IR  -►  IR  ,  equipped  with  the  usual  inner  product  ^x,y^  and 
norm  i|x l|  . 

2 

For  each  T  >  0,  L  j  is  the  set  of  all  periodic  measurable  functions 

2 

x:  IR  -►IR  with  period  T  such  that  the  integral  of  x  over  one  period  is 
finite.  It  is  a  Hilbert  space  with  the  "average  power"  inner  product 

T 

4  Y  I  x(t)'y(t)  dt  .  (14) 

0 


2 

The  norm  on  L  j  is  denoted 


Linearization  of  the  source  admittance  operator  about  a  nominal  input 


waveform  is  conventionally  accomplished  using  the  Frechet  derivative 

[10,  sect.  3.1],  [11,  sect.  2.1c].  Unfortunately,  one  cannot  encompass 

2 

nonlinear  resistors  in  a  framework  that  requires  differentiation  on  L  or 

p 

L  j.  It  is  a  little  known  fact  that  if  f:  IR  IR  and  the  memoryless  operator 

p  p 

F:  x(t)  f(x(t))  maps  L  to  L  ,  then  F  is  not  Frechet  differentiable  unless 

f  =  ax  +  b  [12,  Appendix].  For  this  reason  the  theory  in  this  paper  is  based 
on  the  weaker  Gateaux  derivative  [10,  sec.  4.1],  [11,  sect.  2.1c]. 

Given  L,  L'  as  above,  an  operator  F:  L  L',  and  x,heL,  suppose  there 
exists  an  element  denoted  dF(x,h)  of  L'  such  that 

lim  |l  F(xnh)  -  F(x)  I  -q 

t-O+ji  t  u  . 

Then  5F(x,h)  is  called  the  Gateaux  variation  of  F  at  x  for  the  increment  h 
[10, p. 251].  If  fF(x,h)  exists  for  all  x,heL,  and  if  for  each  xcL  the  map 
h  -  5F(x,h)  is  an  element  of  L(L,L'),  then  F  is  said  to  be  Gateaux 
differentiable  on  L.  In  this  case  the  map  x  6F(x,*)  is  called  the  Gateaux 
derivative  of  F  and  denoted  OF:  L  L(L,L')  [10,  pp. 255-256].  Similarly, 
iF(x,-)  is  denoted  DF(x)eL(L,L'),  and  6F(x,h)  is  denoted  (DF(x))heL'. 

2.2)  Main  Theorem 

Theorem  1  (Maximum  Average  Power  in  Periodic  Steady  State) 

Fix  T  >  0  and  let  in  Fig.  1  be  characterized  by  an  admittance  oper- 
2  2 

ator  F;  Ly  -  Ly,  where  Ij  is  any  linear  subspace  of  L ^  .  Suppose  F  is 
Gateaux  differentiable  on  Lj  and  the  associated  operator  H:  Lj  -►  Ly, 

2.  Thus,  if  v(-)  has  period  T  and  lies  in  Ly  ,  the  steady-state  response  i(-) 
of  cannot  have  subharmonics. 
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given  by 


adj 


H:  V  I-*  F{v)  +  (DF(y))  v 


is  strictly  increasing. 

Then  for  each  i^cH(Lj)  there  is  a  unique  solution  v(ig)eLj  to 
the  network  equation, 

is  =  H(v)  , 

3 

and  the  average  power  absorbed  by  the  load, 

P(v)  A  "  Os-  ''(v),v>.p 

has  a  unique  global  maximum  over  Ly,  which  is  attained  at  v  =  v(ig). 


(15) 


(16) 


(17) 


Corol lary  (Maximum  Total  Energy  for  Transients) 

2 

Let  L  be  a  linear  subspace  of  L  and  substitute  L  for  Lj  in  the 

4 

assumptions  of  Theorem  1.  Then  the  same  conclusions  hold,  but  with 
v(i^)£ L  maximizing  the  total  energy  E(v)  A  <^ig  -  F(v),y^  over  L. 

Note  that  H  is  the  sum  of  the  source  admittance  F  and  the  optimal 
load  admittance  operator 


(18) 


as  stated  less  formally  in  Eq.  (2). 

In  applications  one  might  wish  to  restrict  attention  to  currents  and 
2 

voltages  in  Ly  with  additional  properties  such  as  continuity  or  boundedness 

2 


This  is  the  reason  for  introducing  LyC L y  in  the  formulation  of  Theorem  1 


A  more  explicit,  but  cumbersome,  notation  would  be  ^(v.ij).  Using  it. 
Theorem  1  states  that  Vv,i^£Ly,  P(v,i^)  <  P(v(ig)  ,i^)  if  v  f  v(i^). 

For  the  Corollary,  the  adjoint  is  of  course  taken  with  respect  to  the 
inner  product  on  L^  rather  than  ^  ,  ^y  . 


The  essential  idea  behind  ttie  theorem  is  that  a  solution  v(*)  of  Eq.  (16) 
is  a  stationary  point  of  -"IR,  and  the  monotonicity  assumption  on  H 

guarantees  that  P"  is  strictly  concave.  Details  follow. 


Proof  of  Theorem  1 


Uniqueness  of  the  solution  to  Eq.  (16)  follows  from  the  fact  that  H  is 
strictly  increasing.  By  the  chain  rule  for  the  composition  of  Frechet-  and 
Gateaux-differentiable  functions  [10, p. 253],  F  is  Gateaux  differentiable 
and  for  all  x,h  e  Lj, 

(DP(x))h  =  {i^-  F(x),h)^  -  <(DF(x))h,x)^  = 

^•5  -  F(x)  -  DF(x))^‘^J  x.h)^  =  (19) 

<^■5  -  H(x),h)^  . 

Thus,  if  i.eH(L,) 

5  1  9 

a)  DP(0(i^))  =  OzLiljXj)  , 

b)  given  any  x.yeLy,  the  map  \  ^  P[x+\(y-x)]  is  differentiable  at 
each  \-:]R  ,  and 

c)  ^P[x+-(y-x)]  =  -  H[x+A(y-x)],  (y-x)^y  . 


To  show  that  v(i^)  globally  optimizes  P,  fix  i^=;H(Ly),  let  v  =  v(i^), 
and  choose  any  veLy,  v  v  .  Then 

Rv)  -  P(^)  = 

F[v+'‘ (v-v)]  i  -  F[v+\(v-v)]  I 

'^=1  !a=o 

^  F[v+\(v-v)] }  dX  . 

0 


(20) 


Using  c),  the  integrand  above  is 


-  H[v+\(v-v)].  v-v)^y  = 

(since  i^  =  H(v)) 

^  ^H[v+X(v-v)]  -  H(v),  [v+X(v-v)]  -  [v]^y,  VX>0  , 

and  the  integrand  vanishes  at  \=0.  The  inner  product  above  is  strictly 
positive  for  X?«0  since  H  is  strictly  increasing  by  assumption.  Thus,  the 
integrand  in  Eq.  (20)  is  negative  for  x>0  and  zero  for  X=0,  so  F(v)  <  F(v) 
as  claimed.  | 

The  proof  of  the  Corollary  is  essentially  identical  and  hence  is  omitted. 

2.3)  Iterative  Algorithm  for  Determining  v(0 

Equation  (IS)  shows  that  i^  -  H(v)  is  the  gradient  [10,  p.l96], 

[11,  sect.  2.5],  [13,p.54ff]  of  Pat  v.  This  suggests  that 

we  attempt  to  maximize  ^  by  a  simple  "hill -cl imbing"  algorithm  of  the  form 

Xj^l  =  r(1^-H(Xj))  +  Xj  ^  M(Xj)  (21) 

for  some  r>0.  Note  that  under  the  assumptions  of  Theorem  1,  j_f  {v^}  converges 
to  some  veLj  and  H  is  continuous,  then  i^  =  H(v)  and  v  globally  maximizes  P. 
By  tightening  the  assumptions  a  little  further,  we  can  guarantee  convergence 
for  all  sufficiently  small  positive  r. 

Theorem  2 

Strengthen  the  assumptions  of  Theorem  1  by  supposing  further  that  Lj 

is  closed  and  H  is  uniformly  increasing  and  Lipschitz  continuous  on  Ij. 

2 

Then  for  any  i-eLy,  any  initial  guess  v  eLj,  and  any  r£(0,  25/K  ),  the 
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sequence  generated  by  Eq.  (21)  converges  to  v(ig). 

Note  that  Theorem  2  also  guarantees  existence  of  a  solution  to  Eq.  (16) 
for  all  i.e.,  H(Ly)  =  Lj.  The  proof  is  a  straightfonvard  application 

of  [14]  and  is  similar  to  that  in  [15].  It  is  given  in  detail  in  [2]  and 
will  be  omitted  here. 

Note  that  the  example  in  Sect.  1.2)  satisfies  all  the  assumptions  of 

Theorem  2  except  that  of  global  Lipschitz  continuity,  since  the  derivatives 

g2'(-)  and  93'(')  are  unbounded.  (Because  they  are  bounded  on  every  bounded 

subset  of  IR,  a  more  detailed  argument,  omitted  here,  shows  that  the  solutions 
—  00  2 

obtained  maximize  P  over  LjOLy  ,  which  is  certainly  sufficient  in 
practice. ) 

IV.  Concluding  Remarks 

The  interested  reader  may  wish  to  compare  these  results  with  those  ob¬ 
tained  by  a  describing  function  method  in  [16]. 
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Fiqure  Captions 

Fig.  1  The  optimal  nonlinear  source  admittance,  given  in  Eqs.  (2)  and  (18), 
is  independent  of  the  current  source  waveform  i5(’)  and  extracts 
the  maximum  power  from  the  source  for  every  igC’).  Although 
is  noncausal  in  general,  it  can  be  useful  in  designing  a  realizable 
load  that  maximizes  power  for  some  particular  igC')- 

Fig.  2  Theorem  2  enables  one  to  numerically  determine  the  optimal  output 
voltage  v(-)  for  this  circuit  when  the  resistor  curves  are  as  shown 
in  Fig.  3. 

Fig.  3  The  three  resistor  curves  for  the  circuit  in  Fig.  2  are 
g,^(v)  1  v'v  k=1.2,3,  with  gj(0)  A  0  . 

Fig.  4  One  period  of  the  optim.al  output  voltages  for  the  circuit  in  Fig.  2. 


